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ABSTRACT 

We present a new constraint on the biased galaxy formation picture. Gravitational instability theory 
predicts that the two-point mass density correlation function, £(r), has an inflection point at the sep- 
aration r = r , corresponding to the boundary between the linear and nonlinear regime of clustering, 
£ ~ 1. We show how this feature can be used to constrain the biasing parameter, b 2 = £s( r )/£( r ) on 
scales r ~ r 0l where £ g is the galaxy-galaxy correlation function, allowed to differ from £. We apply our 
method to real data: the £, g {r), estimated from the APM galaxy survey. Our results suggest that the 
APM galaxies trace the mass at separations r ;> 5/i _1 Mpc, where h is the Hubble constant in units 
of 100 km s _1 Mpc _1 . The present results agree with earlier studies, based on comparing higher order 
correlations in the APM with weakly nonlinear perturbation theory. Both approaches constrain the b 
factor to be within 20% of unity. If the existence of the feature we identified in the APM £, g (r) - the 
inflection point near £ 9 = 1 — is confirmed by more accurate surveys, we may have discovered gravity's 
smoking gun: the long awaited "shoulder" in £, predicted by Gott and Rees 25 years ago. 
Subject headings: Cosmology - large-scale structure of universe 
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1. INTRODUCTION 

The concept that galaxies may not be fair tracers of the 
mass distribution was introduced in the early 1980s, in 
part in response to the observation that galaxies of dif- 
ferent morphological types have different spatial distribu- 
tions and that therefore they cannot all trace the mass 
(there are two excellent reviews on the subject: Strauss 
& Willick 1995 and Hamilton 1998). However, there was 
also another reason: to "satisfy the theoretical desire for 
a flat universe" (Davis et al. 1985, p. 391). More precisely, 
biasing was introduced to reconcile the observations with 
the predictions of the Einstein-de Sitter cold dark matter 
(CDM) dominated model. At the time, it seemed that 
just a simple rescaling of the overall clustering amplitude 
by setting £, g (r) = b 2 £,(r), where b « 2 might do the job 
(Davis et al. 1985). However, very soon thereafter, it be- 
came clear that this is not enough: while the unbiased 
(b = 1) £(r) had too large an amplitude at small r, the 
biased model did not have enough large-scale power to 
explain the observed bulk motions (Vittorio et al. 1987, 
Ostriker 1993). A similar conclusion could be drawn form 
comparison of the relative amplitude of clustering on large 
and small scales (eg Maddox et al 1990). The problem 
with the shape of £(r) became explicit when measurements 
of £ g (r) showed that the optically selected galaxies follow 
an almost perfect power law over nearly three orders of 
magnitude in separation. This result disagrees with N- 
body simulations. The standard (£!.,„ = 1) CDM model 
as well as its various modifications, including Q, m < 1 and 



a possible non-zero cosmological constant, fail to match 
the observed power law (see Fig 11-12 in Gaztanaga 1995, 
Jenkins et al. 1998; most of these problems were already 
diagnosed by Davis et al. 1985). Two alternative ways 
out of this impasse were recently discussed by Rees (1999) 
and Peebles (1999). A possible response to the CDM cri- 
sis is to build a model where simple phenomena, like the 
power-law behavior of £ g are much more complicated than 
they seem. In particular, one can explore the possibility 
that the emergence of large scale structure is not driven 
by gravity alone but by "environmental cosmology" a 
complex mixture of gravity, star formation and dissipative 
hydrodynamics (Rees 1999). A phenomenological formal- 
ism, appropriate for this approach was recently proposed 
Dekel & Lahav (1999). There are also some recent analyti- 
cal models, based on halo profiles and halo-halo clustering 
(Seljak 2000, Scoccimarro etal 2001). An obvious alter- 
native to environmental cosmology was recently discussed 
by Peebles (1999), who pointed out that "as Kuhn has 
taught us, complex interpretations of simple phenomena 
have been known to be precursors of paradigm shifts". 
Instead, one can explore a simpler option, that galaxies 
trace the mass distribution, at least for local (low red- 
shift), optically selected galaxies with a broad magnitude 
sampling. This approach rests on the idea that no matter 
how or where galaxies form, they must eventually fall into 
the dominant gravitational wells and therefore trace the 
underlying mass distribution (see Peebles 1980, hereafter 
LSS; Fry 1996). 



The absence of biasing on large (weakly nonlinear) scales 
agrees well with other observational evidence. The mea- 
surements of the two-, three- and four-point connected 
moments of the density field in the APM catalog pro- 
vide support for the hypothesis that galaxies trace the 
mass and also that the large-scale structure we observe 
today grew out of small-amplitude, Gaussian density fluc- 
tuations in an expanding, self-gravitating non-relativistic 
gas. Indeed, the theoretical predictions for the first few 
connected moments, based on this hypothesis (Juszkiewicz 
et al. 1993, Bernardeau 1994) are in good agreement with 
the APM measurements (Gaztahaga 1994, Gaztahaga & 
Frieman 1994, Frieman & Gaztahaga 1999, and references 
therein). The current precision of this higher order cor- 
relation test is 20% and expected to improve with future 
data. The absence of biasing is also suggested by the most 
recent measurements of the mean relative pairwise velocity 
of galaxies (Juszkiewicz et al. 2000). 

In this Letter we propose a new test of the biasing hy- 
pothesis. Our test is based on the behavior of £(r) at the 
£ = 1 nonlinearity boundary. We describe our theoretical 
model in the next section. It is checked against N-body 
simulations in §3. It is then applied to the APM survey 
in §4. Finally, in §5 we discuss our results. 

2. THE INFLECTION POINT 

In the gravitational instability theory, newly forming 
mass chimps are generally expected to collapse before re- 
laxing to virial equilibrium. If this were so, the relative 
mean velocity of particle pairs |i>i 2 (r)| would have to be 
larger than the Hubble velocity Hr to make v 12 (r) + Hr 
negative. As a consequence of the pair-conservation equa- 
tion (LSS), the slope dln£(r-)/dlnr = — 7(7-) must rapidly 
decrease with decreasing separation near the boundary of 
nonlinearity, i.e. when r ~ r . This effect was recognized 
long ago by Gott & Rees (1975). When the expected 
"shoulder" was not found in the observed galaxy-galaxy 
correlation function, Davis & Peebles (1977) introduced 
the so-called previrialization conjecture as a way of reduc- 
ing the size of the jump in 7(7") (the conjecture involves 
non-radial motions within the collapsing clump; see the 
discussion in LSS, §71 and Peebles 1993, pp. 535 - 541; see 
also Villumsen & Davis 1986; Lokas et al. 1996 and Scocci- 
marro & Frieman 1996). Later observational work showed 
a shoulder in £ in several redshift and angular catalogs, 
which was also interpreted as evidence for the boundary 
between linear and nonlinear gravitational clustering (see 
the review by Guzzo 1997 and references therein). 

Quarter a century later the precision of N-body simula- 
tions as well as the quality of the observational data have 
improved dramatically enough to justify a reexamination 
of the problem. The actual shape of the correlation func- 
tion near £ = 1 can be investigated with high resolution 
N-body simulations like those run by the Virgo Consor- 
tium (Jenkins et al. 1998). As shown in Juszkiewicz et 
al. (1999, hereafter JSD), in all four of the Virgo models 
(SCDM, ACDM, OCDM, tCDM) the slope of £(r) exhibits 



a striking feature. Instead of a shoulder, or a simple dis- 
continuity in 7(7"), however, £(r) has an inflection point, 
d 2 £,(r)/dr 2 = which occurs at a uniquely defined separa- 
tion r = r„. At this separation, the logarithmic slope of £ 
reaches a local maximum, rfln^/dlnr = —7*. In all mod- 
els JSD investigated, the inflection point indeed appears 
near the transition £ = 1, as expected by the earlier spec- 
ulations, involving the "shoulder" in £. The separation r* 
is almost identical with the scale of nonlinearity: 
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More precisely, a comparison of Figure 1 in JSD with Fig- 
ure 8 in Jenkins et al. (1998) gives 



\r — r*| < 0.1 r 



(2) 



for all four considered models. Moreover, for all models, 
studied by JSD, the —7 vs. r dependence can be described 
as an S-shaped curve, with a maximum at r = r* ~ r , 
and a minimum at a smaller separation. The depth of 
the minimum in —7(7") increases with increasing normal- 
ization parameter, erg - the final linear rms mass density 
contrast, measured in spheres of a radius of 8 /i _1 Mpc (see 
Fig. 8 in JSD). If the relation P) is indeed a general prop- 
erty of gravitational clustering, it can be used as a test 
of biasing as follows. Suppose the biasing factor is signif- 
icantly greater than unity: b 3> 1. Then £ ff 3> £ and the 
relation (pi) will break down. For power-law galaxy cor- 
relation function, £ 9 (r) = {r og /r)' 1 = b 2 £(r), and instead 
of equation (|l|) we will have r* w r og b^ 2 ^. Since the 
observed slope is 7 sw 1.8, for b = 2, the shoulder in the 
correlation function should appear at a separation smaller 
than a half of the r og parameter! The technique we pro- 
pose is unable to constrain more baroque biasing models 
with a large number of free parameters. However, the pre- 
dictive power of such models is questionable and one may 
ask: are they falsifiable at all? 

3. APM-LIKE N-BODY SIMULATIONS 

In this section, we compare the JSD analytic and N- 
body results with a different set of P 3 M simulations. In- 
stead of the family of CDM models, considered by JSD, 
we use APM-like initial conditions, which have Gaussian 
initial conditions with an initial power spectrum designed 
to evolve into a final spectrum, matching the APM mea- 
surements, under the additional assumption of no bias and 
fl m = 1. We use simulations with identical APM-like spec- 
trum, with A = and two different values of the density 
parameter: f2 m = 1 and 0.3. The box size is 600 h^ 1 
Mpc (or 300 h- 1 Mpc) with 200 3 (or 100 3 ) dark matter 
particles with er 8 — 0.85 (for more details see Baugh & 
Gaztahaga 1996). 

The evolved, nonlinear correlation functions, measured 
from simulations are shown in Figure u\ (top left panel) . 
The full squares correspond to the Q m = 1 model, while 
the open squares represent O m = 0.3. For comparison, we 
show the linear correlation function (dashed line). Nonlin- 
ear effects are more pronounced in the low density model. 
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Fig. 1. — Top left panel: linear f(r) (dashed line) and the measured nonlinear £(r), obtained from the APM-like simulations with Q m = 0.3 
(open squares) and fi m = 1.0 (full squares). Bottom left panel: Corresponding logarithmic slope, — "f(r) = dln^/dlnr for each of the three 
curves from the top panel. The right panels show similar results for the spatial §(r) estimated for APM galaxies (symbols with errorbars), 
compared to the same linear theory APM-like model (dashed line). The vertical dotted lines shows the scale r , defined by the condition 
£( r o) = 1 (top) and the scale r*, at which the nonlinear 7(r) curve crosses the linear one (bottom). 



Note however, that although the correlation functions dif- 
fer significantly in amplitude at separations r < 2ft,~ 1 Mpc, 
their slopes j(r) are almost indistinguishable. 

The particle resolution (the Nyquist wavelength ex 
TV -1 / 3 ) of the simulations used here is significantly lower 
than the resolution of Virgo simulations, and the noise 
in the measured £(r) is further amplified by differentiat- 
ing over r. As a result, determining the position of the 
inflection point r* directly from the j(r) curve alone is 
more difficult. A noise-resistant, alternative definition of 
r* , suggested by JSD is to identify r* with the separation 
at which the linear 7(7"), derived from initial conditions, 
crosses the nonlinear 7(7-) curve, measured from the simu- 
lations. This approach is indeed effective for our APM-like 
simulation (see Fig.f ), but it can not be applied to the true 
APM data since the we know only the nonlinear 7(7-). The 
observed — 7(7") curve does resemble the S-shaped —7(7-) 
from the VIRGO simulations; the main difference is that 
the peak as well as the trough are broad and fuzzy rather 
than narrow and sharp as in the simulations. To deal with 
this problem, we used the following prescription. For a 
discrete set of measurements 7*i, 7 (7^), i — 1 , 2, . . . starting 
with some separation 7*1 rightward of the peak, and moving 
to smaller separations, we compare consecutive values of 
—7 (ft) and — 7(7-4+1). We then identify r* with the largest 
Tj, for which —7(7^) drops below its maximum value for the 
next few points (this is to avoid picking up local maxima 
due to fluctuations, see Fig.l). Obviously, the above pre- 
scription would also pick r* in high resolution simulations. 
When applied to the APM-like simulations, this method 
gives r* ~ 4 — 5 h" 1 Mpc. The same simulations also give 
r ~ 5/i _1 Mpc, in excellent agreement with equation (El), 



which we will consider as a measure of a systematic er- 
ror, introduced by the theoretical model we use. From 
these simulations we conclude that equality between r* 
and r can probably be considered as a generic outcome of 
gravitational dynamical evolution in a model where galax- 
ies trace the mass and the initial slope, dln^/dlnr, is 
a smooth decreasing function of the separation r (as ex- 
pected in hierarchical clustering models, see e.g. LSS). 

4. COMPARISON WITH OBSERVATIONS 

The measurements of £ g (r), and 7(7*) = — dln£ ff /dlnr, 
obtained from the angular correlations of galaxy pairs in 
the APM catalog (Baugh 1996), are plotted in Figure 
\\\ Errorbars correspond to the dispersion in 7(7") from 
4 APM strips of about 30x60 degrees each, eg 200x400 
h~ 1 Mpc at the mean depth of the APM. At scales of 
r ~ 5/i _1 Mpc -C 200 /i Mpc most of the covariance be- 
tween consecutive bins in £ 9 is due to large scale density 
fluctuations that shift the mean density between strips. 
This introduces covariance in the amplitude of £ g (r) (the 
whole £ g curve is shifted vertically from strip to strip, see 
eg Fig. 2 in Baugh 1996), but does not affect much its shape 
7(7-). Thus, we consider our errors in 7(7") as independent. 
The top left panel shows the two-point function (points 
with error bars), and the linear theory curve, described in 
§3 (dashed line). The intersection of the two perpendicular 



dotted lines marks the point (£ g 



(l,r og ). The bot- 



tom panel of Figure G] shows the APM 7(7") as a function of 
the pair separation r. Note the remarkable similarity be- 
tween the empirical data and the characteristic peak in the 
—7(7*) found in the simulations (compare the left and right 
bottom panels in Fig.l; see also Fig.l in JSD). The inter- 



section of the two mutually perpendicular, dotted lines in 
the bottom panel of Figure fil marks the result of apply- 
ing our prescription for estimating r* to the APM data. 
The slope —7 drops down from its maximum value at the 
separation r ~ 5ft. _1 Mpc, and to first approximation this 
scale could be identified with r* . Taking into account the 
error bars in Figure [y we obtain r* ~ (5.5 ± 1.5) h^ 1 Mpc, 
r og ~ (4.5 ± 0.5) h~ 1 Mpc and 7* ~ —1.4. From these 
measurements one can estimate b in the linear bias model 



b 2 = (r og /r*y*. 



(3) 



This expression for b 2 is also exactly valid for a non- 
linear scale-dependent bias at b(r ) as far as r* ~ rn: 
b 2 {r ) = £ ff (r )/£(r ) = {r og /r )^°) = (r og /n)^. This 
is true even if the slopes of the galaxy and matter correla- 
tions are different at r — rg. We find 



6(r )~ 1.15 ±0.23 (±0.11) 



(4) 



at one-sigma level in the errors. The error in parenthesis 
corresponds to the systematic uncertainty in Eq.g. 

5. DISCUSSION 

Recently, Hamilton and Tegmark (2000) found no ev- 
idence of an inflection at the linear-nonlinear transition 
scale in £ ff , estimated from the PSCz survey. The origin 
of this difference with our APM results is not clear, since 
a recent study of the three-point correlations in the PSCz 
catalog (Feldman et al. 2001) leaves little room for biasing, 
providing constraints on b similar to those obtained here 
and to those obtained earlier from the measurements of 
three-point correlations in the APM catalog. If the APM 
inflection as well as its absence in the PSCz are both real 



phenomena, it then means that APM galaxies are less bi- 
ased than PSCz (IRAS) galaxies. Apart from constraining 
the linear bias, Feldman et al. have also measured a small 
but statistically significant second-order biasing parame- 
ter, consistent with the observation that infrared-selected 
(IRAS) galaxies avoid high density cores of clusters. We 
plan a systematic study of this effect in near future. 

We are impressed how well the shape of 7(7*) in the APM 
observations resembles gravity's "shoulder" . This feature 
is a robust result from the APM catalogue and can be seen 
directly in the angular 2-point function and in the recov- 
ered shape of the power spectrum (Baugh & Efstathiou 
1993). Numerical simulations (Gaztanaga & Baugh 1996) 
show that this is not an artifact of the de-projection. It is 
difficult to imagine how such an agreement could happen 
by a mere coincidence, which would have to be the case if 
£ g is unrelated to £. Our results are by no means final, they 
are also less rigorous than one could wish because we are 
limited by the accuracy of the present observational data. 
New generation of catalogs promise an improvement on 
this front in the near future (for an excellent collection of 
reports on the state of the art in this field, see Colombi et 
al. 1998). 
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